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In  our  previous  article  [D.  H.  Huang,  A.  Singh,  and  D.  A.  Cardimona,  J.  Appl.  Phys.  87,  2427 
(2000)],  we  explained  the  experimentally  observed  zero-bias  residual  tunneling  current  [A.  Singh 
and  D.  A.  Cardimona,  Opt.  Eng.  38,  1424  (1999)]  in  quantum-well  photodetectors  biased  by  an  ac 
voltage.  In  this  article,  we  extend  our  theory  to  include  the  photoemission  current  and  reproduce  our 
recent  findings  on  the  dynamical  drop  of  photoresponsivity  7£p h(t)  from  its  static  value  7£ph  in 
quantum-well  photodetectors  as  a  function  of  the  chopping  frequency  of  the  incident  optical  flux.  In 
this  theory,  we  derive  a  dynamical  equation  for  a  nonadiabatic  space-charge  field  S^t)  in  the 
presence  of  an  applied  electric  field  £b(t)  and  an  incident  optical  flux  d>op(f).  From  it,  a 
compensation  of  the  charge  fluctuations  in  quantum  wells  is  predicted  as  a  result  of  dual  tunneling 
and  photon-assisted  escaping  paths.  We  also  find  a  suppression  of  the  nonadiabatic  deviation  of 
7£ph(0  from  71  ph  due  to  a  charge-depletion  effect  in  the  quantum  wells.  ©  2001  American  Institute 
of  Physics.  [DOI:  10.1063/1.1351867] 


I.  INTRODUCTION 

In  two  recent  articles,1*2  we  found  a  residual  tunneling 
current  in  multiple  quantum  wells  when  an  ac  bias  voltage 
sweeps  through  zero.  A  circuit  model1  including  an  impor¬ 
tant  tunneling  resistance  in  series  with  a  quantum-well  ca¬ 
pacitance  was  devised  to  explain  this  phenomenon,1  and  a 
satisfactory  numerical  simulation  was  obtained  by  using  this 
phenomenological  model.  It  indicates  that  a  physical  process 
with  a  very  large  time  constant  is  involved  in  transport 
through  multiple  quantum  wells  (MQWs).  The  microscopic 
origin  of  this  observation  was  explored  thereafter,2  and  a 
current  instability  and  hysteresis,  as  well  as  a  current  “arch” 
and  “ripple,”  were  predicted  and  confirmed  experimentally. 

It  is  well  known  that  resonant  electron  tunneling  in 
MQWs  can  occur  only  when  the  barrier  between  adjacent 
quantum  wells  is  thin.  If  the  banier  is  very  thick,  the  phase 
of  the  wave  function  will  be  completely  lost  as  an  electron 
tunnels  from  one  well  to  another.  As  a  result,  only  sequential 
electron  tunneling  exists  for  thick  barriers.  If  a  dc  electric 
field  £b  is  applied  to  the  system,  electrons  in  quantum  wells 
simply  respond  to  it  through  an  adiabatic  tunneling  current 
fl^b]  which  is  a  nonlinear  function  of  £b  due  to  sequential 
electron  tunneling.  We  have  found  that  when  a  time- 
dependent  electric  field  £b(t )  is  applied  to  the  system,  it 
induces  a  fluctuation  in  the  charge  density  inside  the  quan¬ 
tum  wells  around  the  equilibrium  value  n2 d*  This  gives  rise 
to  a  nonadiabatic  space-charge  field  £m(t)  which  modifies 
the  adiabatic  tunneling  current  It[£b(t)]  by  adding  a  nona¬ 
diabatic  correction  A/f(r).  Under  this  situation,  A It(t)  re¬ 
mains  in-phase  with  £na(t ),  and  the  dynamics  of  £na(f)  are 
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determined  by  the  source  term  d£b(t)ldt  and  the  usual 
quantum-well  char ging/ discharging  process.  If  £m(t)  is 
positive,  which  shifts  the  Fermi  energy  down,  the  quantum 
well  is  discharged  with  its  transient  charge  density  lower 
than  n1D .  The  quantum  well  can  also  be  charged  when  £na(r) 
becomes  negative. 

When  one  uses  a  quantum- well  photodetector  to  look  for 
a  distant  target  buried  in  cold  outer  space  (—4  K),  the  device 
temperature  must  be  kept  very  low  (Te~ 40  K)  in  order  to 
minimize  the  noise  and  enhance  the  signal-to-noise  ratio.  In 
addition,  when  the  target  is  moving,  a  multiple  sampling 
process  is  required  to  detect  the  target  motion  and  reduce  the 
noise  by  turning  the  shutter  of  a  photodetector  on  and  off. 
However,  a  drop  of  the  photoresponsivity  of  the  device  was 
found3  when  the  shutter  frequency  exceeded  a  threshold 
value.  This  threshold  frequency  depended  on  the  device 
temperature,  the  external  bias  voltage  and  the  incident  opti¬ 
cal  flux.  The  effect  of  a  shutter  can  be  simulated  by  a 
chopped  incident  optical  flux  <J>op(f).  When  the  quantum- 
well  photodetector  is  exposed  to  <J>op(0,  the  charge  density 
in  the  quantum  wells  again  fluctuates  around  «2d*  As  ex¬ 
plained  in  the  tunneling  case  earlier,  a  nonadiabatic  space  - 
charge  field  £m(t)  will  be  induced  in  the  system.  Here,  the 
dynamics  of  <£^(0  are  determined  by  the  source  term 
d®op(t)/dt  and  the  usual  quantum- well  charging/dis¬ 
charging  process  with  its  decay-time  depending  on  <*>oP(0. 
Moreover,  £m(f)  not  only  modifies  the  adiabatic  photoemis¬ 
sion  current  Ie[£b& op(0]  by  subtracting  an  out-of-phase 
correction  A Ie(t)  relative  to  4>op(f),  but  also  modifies  the 
tunneling  current  by  adding  an  in-phase  correction  A I  ft). 

When  the  charge  fluctuates  in  the  quantum  wells,  the 
photoresponsivity  of  the  detectors  gains  a  nonadiabatic  de¬ 
viation  from  the  adiabatic  value  occurring  when  no  charge 
fluctuation  (CF)  is  present.  This  will  cause  a  deformation  in 
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the  detected  images.  We  know  that  the  quantum-well  CF  can 
be  individually  controlled  by  either  d<J>op(r)/dr  for  a  dc  elec¬ 
tric  field  or  d£b{t)ldt  with  no  incident  photons.  When  both  a 
time-dependent  electric  field  £b(t)  and  a  time-dependent  in¬ 
cident  optical  flux  5>op(0  are  applied  to  the  system,  the  CF 
in  the  quantum  wells  will  be  determined  by  d£b{t)ldt  and 
dQ>op(t)/dt  simultaneously.  In  this  case,  both  the  tunneling 
and  photon-assisted  escape  channels  are  open  to  electron 
transport.  In  order  to  minimize  the  image  deformation,  we 
need  to  maintain  the  phase  of  £b(t)  opposite  to  that  of 
*op(0.  Under  this  condition,  the  CF  from  these  two  sources 
will  compensate  each  other  in  the  quantum  wells,  and  both 
A It(t)  and  A Ie(t)  can  be  greatly  reduced.  As  a  result,  the 
photoresponsivity  of  the  device  will  approach  its  adiabatic 
value  and  the  system  will  behave  close  to  an  ideal  adiabatic 
one. 

The  organization  of  this  article  is  as  follows.  In  Sec.  n, 
we  present  our  model  beyond  the  adiabatic  limit  by  deriving 
a  general  nonlinear  dynamical  equation  for  the  nonadiabatic 
space-charge  field  in  the  presence  of  both  a  time-dependent 
electric  field  and  a  time-dependent  incident  optical  flux.  For 
the  special  case  with  a  dc  electric  field  and  a  time-dependent 
incident  optical  flux,  the  experimentally  observed  drop  of  the 
photoresponsivity  as  a  function  of  the  chopping  frequency  of 
the  optical  flux  is  reproduced.  Numerical  results  and  discus¬ 
sions  are  given  in  Sec.  EH  for  the  nonadiabatic  space-charge 
field,  total  nonadiabatic  photoemission  current,  and  dynami¬ 
cal  photoresponsivity  when  both  the  electric  field  and  the 
incident  optical  flux  are  time  dependent.  The  article  is  finally 
concluded  in  Sec.  IV. 

II.  MODEL  AND  THEORY 

In  this  section,  we  will  first  study  the  CF  resulting  from 
tunneling  transport  in  the  presence  of  a  time-dependent  elec¬ 
tric  field  £b{t)  and  no  incident  photons.  Next,  the  CF  result¬ 
ing  from  photon-assisted  escape  will  be  explored  when  the 
only  time  dependence  arises  from  an  incident  optical  flux 
<I>Op(0.  Finally,  the  compensation  of  the  CF  in  quantum 
wells  will  be  investigated  when  both  £b(t)  and  <frop(0  are 
present. 

A.  Tunneling 

In  order  to  introduce  notations  and  make  a  comparison 
between  the  CFs  resulting  from  either  the  tunneling  or 
photon-assisted  escape,  we  begin  by  deriving  some  of  the 
equations  in  our  previous  article.2 

Let  us  first  consider  the  tunneling  transport  of  electrons 
in  a  MQW  system  under  a  bias  field  £b(t).  We  find  that  the 
tunneling  current  depends  not  only  on  £b(t)  which  produces 
a  sequential  tunneling  current,  but  also  on  d£b{t)fdt?  In  the 
nonadiabatic  limit,  the  charge  density  in  each  quantum  well 
fluctuates  around  «2d-  It  results  in  a  nonadiabatic  space- 
charge  field  £na(  t)  which  can  be  either  positive  or  negative 
when  the  charge  density  in  the  quantum  wells  is  lower  or 

higher  than /22d- 

To  derive  the  dynamical  equation  for  £na(t),  we  use  Le¬ 
vine’s  sequential  electron  tunneling  model4  to  write  down 
the  adiabatic  tunneling  current  It[£b(t)]  under  the  influence 
of  £b(t)  in  MQWs  as 


It[Sb(t)]=eSvd[£b(t)]ne ^£b{t)Jel 


(1) 


where  Te  is  the  electron  temperature  (or  the  device  tempera¬ 
ture  under  thermal  balance),  S  is  the  sample  cross-sectional 
area,  vd[£b(t)~\  is  the  electron  drift  velocity  related  to  £b(t) 
by  a  saturation  model,  and  ncf££b(t),  Te]  is  the  effective 
three-dimensional  tunneling  electron  density.  In  the  satura¬ 
tion  model,  vJi£b(t)]  is  given  by 


vism= 


'  Vs 

£b(t) 

Vi+[4(0/^]2 

Ss  J’ 

(2) 


where  vs  and  £s  are  the  saturation  velocity  and  field,  respec¬ 
tively.  Moreover,  we  have  defined  in  Eq.  (1): 


nA^b(t)Je]  = 


m* 


7 rh2L 


W! 


r 


dET[E,£b(t)] 


x|/o 

fo 


£+£i-/4c(rj 


kBTe 

£+ E  i  +  eLB£b{t)  —  fic{T e) 


kBTe 


(3) 


where  fo(X)  is  the  Fermi-Dirac  distribution  function  and 
fJLc{Te)  is  the  chemical  potential  of  electrons  in  each  quan¬ 
tum  well.  When  £b(t)  is  applied,  the  electrons  in  the  quan¬ 
tum  wells  will  produce  a  steady-state  current  flowing  in  the 
system.  Although  the  ground-state  wave  function  of  the  elec¬ 
trons  inside  the  quantum  wells  can  be  modified  by  £b(t ),  the 
electron  density  in  each  well  will  not  change.  As  a  result, 
becomes  independent  of  £b(t)  and 
thus  is  simply  denoted  by  Ex-juLc(Te).  In  Eq.  (3),  m*  is  the 
effective  mass  of  electrons,  Lw  is  the  width  of  the  quantum 
well,  and  LB  is  the  thickness  of  the  barrier  between  adjacent 
quantum  wells.  Ex  is  the  ground-state  energy  evaluated  at 
5^(0  — 0,  and  l[E,£b(t)]  is  the  transmission  coefficient  of 
electrons  with  incident  energy  E  through  a  barrier  biased  by 
£b{t).  The  difference  of  the  Fermi-Dirac  distribution  func¬ 
tions  in  Eq.  (3)  comes  from  the  requirement  of  an  occupied 
initial  state  in  one  well  and  an  unoccupied  final  state  in  ad¬ 
jacent  well  for  the  sequential  tunneling  process. 

When  d£b{t)!dt± 0,  there  exists  a  surge  tunneling  cur¬ 
rent  Ist(t)  flowing  out  of  the  quantum  wells  in  addition  to 
Iie>m2  working  in  the  nonadiabatic  limit,  Is, (t)  is 
given  by 


Ist(t)  =  —eS lira  t-  «2d 
Af— >0  al  I 


/  m* 


*  \  r+co 

F  I  dE 


X/o 


-~fo 


\  7r^2, 

E+Ex  +  eLBAt  d£b(t)/dt-fjic(Te) 


kBTe 


Ex-fiLc(Te) 


kBTe 


m*e2SLB\  d£b{t) 

^  7rh2  J  dt 


(4) 


The  existence  of  Ist(t)<*d£b{t)ldt  causes  the  imbalance  of 
the  tunneling  current  flowing  into  and  out  of  a  quantum  well, 
which  drives  the  charge  density  away  from  «2d-  The  charge 
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fluctuation  A  Q(t)  in  the  quantum  wells  induces  a  nonadia- 
batic  space-charge  field  £na(0*  This  gives  rise  to  a  nonadia- 
batic  correction  to  the  tunneling  current 

A It{t)=eS{vd[£b{t)+£J<t)]  nefl£b(t)+£Jt),Te] 

-vl£b{t)-]ntfl£b{t),Te1}.  (5) 

In  terms  of  £na(t ),  AQ(t)  in  each  quantum  well  can  be  ex¬ 
pressed  as 


AQ{t)  =  eS\  n2D- 


X 


r 

Jo 


dEf0 


E+El  +  eLB£m(t)-fic(Te ) 


kBTe 


(6) 


Consequently,  the  quantum-mechanical  continuity  equation 
Alt(t)+Ist(t)  +  dAQ(t)/dt=0  leads  us  to  the  following 
dynamical  equation  for  £m(t): 

d£Jt) 

cQw[4,(*)]-^ 

d£b(t)  I  eS\ 

=  CQw[°]  + 

"^"^na(0  V  ^£b{t)^n  Qf^£b{t)  yT y  (7) 

where  CQW[£na(/)]  is  the  dynamical  quantum-well 
capacitance5  given  by 


CQw[£naW]-/o 


E1 +eLB£ja(t)  —fic(Te) 

/  m*e2S\ 

kBTe  .J 

trh2 

(8) 

which  is  different  from  the  QW  capacitance  CQW[0].  If 
and  /j.c(Te)/eLB  for  a  slowly  varying  £b(t), 
we  can  expand  Eq.  (7)  to  first  order  in  £m(t)  and  arrive  at  a 
linear  approximation 

£JLt) 

(9) 


d£m(t)  d£b(t) 


dt  dt  /a4M]cQw[or 

where  Rj_£b{t)']  is  the  differential  tunneling ’  resistance, 
given  by 

w,w]=y  d£~b{vi£b(t^  re]} . 

(10) 

Here,  Te  is  kept  constant,  and  its  dependence  is  not  explicitly 
written  out.  Equation  (9)  is  a  circuit  equation  with  respect  to 
£na(t)  in  the  presence  of  a  source  term  d£b(i)ldty  in  which 
^tr[^(0]^Qw[0]  plays  the  role  of  a  charging! discharging 
time  constant.  This  uncovers  the  microscopic  origin  of  our 
previous  circuit  model1  and  explicitly  relates  the  charging! 
discharging  time  constant  to  the  change  of  microscopic  tun¬ 
neling  current  in  MQWs. 

Using  Eq.  (9),  we  have  predicted  a  current  hysteresis  and 
a  current  arch  for  a  sinusoidal  £b(t),  as  well  as  a  current 
ripple  for  a  step-like  £b(t).  A  current  instability  is  also  found 
by  using  Eq.  (7).  All  of  these  predictions  have  been  con¬ 
firmed  by  our  previous  experiments.1,2 


B.  Photon-assisted  escape 

In  the  presence  of  incident  photons,  electrons  in  the 
ground  state  can  transit  to  the  upper  excited  state  by  absorb¬ 
ing  photons.  From  this  excited  state,  they  can  easily  tunnel 
out  to  the  continuum  states  above  the  barrier  with  help  from 
a  dc  electric  field  £b .  Consequently,  a  photon-assisted  escape 
channel  is  opened  for  electrons  to  get  out  of  quantum  wells 
in  addition  to  the  previous  tunneling  channel. 

In  the  adiabatic  limit,  by  using  Levine’s  electron  photo¬ 
emission  model4,  we  can  write  the  adiabatic  photo-emission 
current  as 

Ie\-£b  > >®op(*)]  =  ^  £ ^]rt2D^)op(0  >  (H) 

where  co  is  the  frequency  of  the  incident  light,  Pe[£b]  is  the 
escape  probability  of  electrons  from  the  upper  excited  state 
to  the  continuum  states  above  the  barriers,  o-0p|>,  £b]  is  the 
optical  cross  section  which  is  related  to  the  absorption  coef¬ 
ficient  by  crop[ct>,  £&]  =  /3abs[<o,  £b]Lw/n2D  in  the  limit  of 
£abs[w>  £b\Lw<\,  and  $op(t)  is  the  incident  optical  flux. 
For  the  escape  probability,  we  use  the  following  empirical 
formula4 


I +A0  expl  - 


^es 


-1 


(12) 


where  A  0  is  the  zero-field  escape  time  ratio  and  £^  is  the 
effective  bamer  lowering  field.  A  more  accurate  escape 
probability  can  be  calculated  by  using  the  time  evolution 
method.6  However,  the  use  of  Eq.  (12)  is  adequate  for  eluci¬ 
dating  the  basic  physics  for  the  compensation  of  CF  in  quan¬ 
tum  wells.  For  the  absorption  coefficient,  we  have7 

£abs!>>  4]=  ^ 


nr[co,  4]\c( 
xll+Pph(h<»/kBTe)]Im  aL[(o,£b]  , 

(13) 

where  pph(2f)  is  the  Bose-Einstein  distribution  function  for 
photons.  In  Eq.  (13),  eb  is  the  relative  dielectric  constant  of 
quantum  wells, 

nr[<*>,  £*]  =  [£+  5  Re  aL[<o,  £b] 

+  |  VO+Re  aL[aj,  £b])2 

+  |(ImaI[w,4])2]1'2  (14) 

is  the  dynamical  refractive  index  function,  and  the  Lorentz 
ratio  is  given  by 

.2 

ll(£i(z)|z|£o(z)>l2 


«lO,  £b\-- 


n  2D* 


€0  €bLW 

1 


1 


(h<o-h£ll0[£b])  +  iy  (h w  +  ^ft10[^])  +  i y\  9 

(15) 

where  hO,i0[£b]  is  the  energy  separation  between  the  ground 
and  upper  excited  states,  y  is  the  homogeneous  energy-level 
broadening,  and  \<£i(z)\ez\£Q(z)>\2  is  the  square  of  the 
transition  dipole  moment  between  the  ground  state  £0(z)  and 
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excited  state  ^  (z).  Here,  g0(z)  and  £i(z)  depend  on  £b  due 
to  the  Stark  effect  In  Eq.  (15),  the  Coulomb  renormalization 
of  the  electron  energy  levels  can  be  included  using  the  self- 
consistent  Hartree-Fock  calculation.7  The  many-body  depo¬ 
larization  effect7  has  been  neglected.  It  will  shift  the  absorp¬ 
tion  peak  slightly  due  to  the  screening  of  the  Coulomb 
interaction  between  electrons. 

Because  4>op(0  varies  with  time,  it  induces  a  CF  in  the 
quantum  wells  in  the  nonadiabatic  limit,  once  again  giving 
rise  to  a  A Q(t)  and  a  A It(t)  as  given  in  Eqs.  (5)  and  (6).  In 
this  case,  however,  we  have  an  additional  nonadiabatic  cor¬ 
rection  to  the  photoemission  current  flowing  out  of  the  quan¬ 
tum  wells,  given  by 

AIe(t)=eSPe[£b]crop[<o,£b]  <I>op(r) 


X/o 


£  +  ffl  +  g£g£na(Q-,Mr<>) 

kBTe 


(16) 


Here,  we  have  neglected  the  secondary  corrections  due  to 
£na (0  t0  the  escape  probability  and  optical  cross  section. 


Moreover,  because  d<t>op(t)Idt=fc  0,  there  exists  a  surge  es¬ 
caping  current  Ise(t)  flowing  out  of  quantum  wells 

d<&  ft) 

/*(0  =  -eS^Kp[w,^]n2D  T,  — jp,  (17) 

where  r,  is  the  lifetime  of  the  excited  electrons.  The  exis¬ 
tence  of  Ise(t)<xd®op(t)/dt  induces  the  imbalance  between 
the  emission  and  capture  currents  flowing  out  of  and  into  a 
quantum  well,  which  deviates  the  charge  density  away  from 
n2D.  The  capture  probability  in  Eq.  (17)  is  given  by  the 
empirical  formula4 


PcVSt\= 


1 4  exp 


(18) 


with  Bm  and  £cp  being  the  high-field  capture  coefficient  and 
effective  well  capturing  field,  respectively.  Capture  probabil¬ 
ity  can  be  more  accurately  calculated  by  adopting  the  Fermi 
golden  rule  method.8  Applying  the  quantum-mechanical  con¬ 
tinuity  equation  A/,(f).+  A/e(t)  +Ise(t)  +  dAQ(t)/dt =  0, 
with  A It(t)  and  A Q(t)  given  by  Eqs.  (5)  and  (6),  leads  to  the 
following  dynamical  equation  for  £na(t): 


J 


d£n  Jit)  leS' 


d®  opM 


CQW[£na(0]  1  JZ  PeVblcTo p[«>,  ^]«2D  T, 


dt 


e£ 

Lb 


+  £na(t) >  Te\  v (K.^b\  nef^L^b  >  P ei} 


j  m*eS\ 

\7 rh^L* j  *oP« 


y-\l‘dE\u 


E+E1-p,c(Tey 

_  f 

'E+E,  +  eLB£J,t)-ixc{Te)V\ 

kBTe  \ 

JO 

kBTe  jj 

(19) 


In  the  limit  of  \£a^t)\<\£b\  and  /xc(Te)/eLB  for  a  slowly 
varying  <J>op(0>  we  can  expand  Eq.  (19)  to  first  order  in 
£„,,(/)  and  arrive  at  a  linear  equation 

+ rJa,  *<,(*)]}  ’ 
(20) 

where  the  optoresistance  and  phototraverse  time  are  defined, 
respectively,  by 

Rop[£b ,  <bop(r)]  =P*t£d°' opt"’  £bl  *Jt)  CQW[0]  , 

(21) 

ropl£b]  =  Tt  Pe[£b\ crop[  a),  4]«2d  ■  (22) 


In  Eq.  (21),  l//?op[£fr,  <I>op(f)]  describes  how  easy  one 
electron  can  transit  from  the  ground  state  to  the  upper 
excited  state  and  then  escape  out  to  the  continuum  state 
above  the  barrier.  By  comparing  Eq.  (20)  with  Eq.  (9),  we 
get  the  effective  resistance  \!ReS^£h,  <bop(0]  =  l/^tr[4] 
+  l/Rop[£b  ,<I>op(/)]  in  the  circuit  equation  for  the  dual  tun¬ 
neling  and  photon-assisted  escape  channels,  the  source  term 
(e'5'7‘op[^]/f'sC’QW[0])dC>op(f)/rfr,  and  the  charging 1 
discharging  time  constant  R^£b  ,<£>op(0]C’qW[0]. 


C.  Compensation  of  charge  fluctuations 

In  the  presence  of  both  £b(t)  and  <l>op(t),  from  the 
quantum-mechanical  continuity  equation  AI,(t)  +  Is,  (t) 
+  Ale(t)+Ise(t)+dAQ(t)/dt=0  we  get  the  dynamical 
equation  for  £m(t)  by  combining  Eqs.  (7)  and  (19): 
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d£JJ)_r,  _ft_  dSb(t)  te&rJAMl]  d<t>op(t) 


Qjwt^iiaW]  ~J~t  —  Cqw[0]  H  ]T 


J^]{vA£b(t)+£w(t)l  nt{l£h(t)+£n,(t),Te]-vd[£b(t)]  nef££b(t),  Te}} 


TTh2LB j  (^op[£fcW>  ^,op(0]CQw[0]J 

■+“  J  lE+Et-fj-ATe)}  „  \E+Ei+eLB£n3(t)~iic(Te) 


XJo  "W 


In  the  limit  of  slowly  varying  £b(t )  and  <£op (t),  by  combin¬ 
ing  Eqs.  (9)  and  (20)  we  are  led  to  the  linear  approximation 
ofEq.  (23): 

CqwCO]  ^ 

rrtl  d£b(t)  t  ( e5Top[fi(f)])  d<S>op(t) 

-cQW[0]  -#-+[  -j-  |  — 

1  1  ) 

~£Jt)  j  wwl + *<»(*)])’  (24) 

The  individual  source  terms  d£b(t)Idt  and 
{e«S7-op[4(z)]/^BCQW[0]}J<I)op(f)/^  in  Eq.  (24)  can  be 
either  in-phase  or  out-of-phase  with  each  other.  When  they 
are  in-phase,  the  effects  of  the  CF  from  these  two  channels 
become  constructive.  When  they  are  out-of-phase,  the  effects 
are  destructive.  This  gives  rise  to  a  compensation  of  the  CF 
in  quantum  wells  from  dual  tunneling  and  photon-assisted 
escape  channels. 

Since  the  general  nonlinear  Eqs.  (7)  and  (19)  and  their 
linear  approximations  in  Eqs.  (9)  and  (20)  are  all  confirmed 
by  our  experiments,  the  validity  of  Eq.  (23)  by  combining 
Eqs.  (7)  and  (19)  or  its  linear  approximation  in  Eq.  (24)  by 
combing  Eqs.  (9)  and  (20)  is  guaranteed.  Therefore,  we  be¬ 
lieve  that  the  prediction  made  from  Eqs.  (23)  and  (24)  in  our 
numerical  results  later  for  the  compensation  of  the  CF  in 
quantum  wells  should  be  observable  in  a  future  experiment. 


D.  Photoresponsivity 

The  dynamical  photoresponsivity  of  the  system  in  the 
presence  of  £b(t)  and  <&op(f)  is  defined  by  the  photoemission 
current  with  transient  charge  density  at  time  t  in  each  quan¬ 
tum  well4  and  is  given  by 

h<oPc[£b(t)]S<t>0p(t)  '  {  > 

In  the  special  case  of  d£b{t)fdt— 0  and  a  slowly  varying 
<&op(0,  we  find  from  Eq.  (25) 


nph(t)  M1  \ePel£b]<rJL<o,£b]n  2D4U05 


^op >  *CO<0] 


where  the  static  photoresponsivity  is 


7?  0  = 

Kph 


ePe[£fc]cropl>,  £fe]»2D 

hcoPc[£b] 


By  assuming  an  optical  flux  <I>  op(0  =  <J>0 + 3>mexpO'fV), 
where  ft  c  and  4>m  are  the  frequency  and  modulation  ampli¬ 
tude  of  <t>op(r)  and  <D0  is  the  background  optical  flux,  we 
take  the  Fourier  transform  of  Eq.  (20)  in  the  limit  of  weak 
modulation  3>m/<I>0<^l  and  get  the  leading  term  of  <fna(ftF); 


eST^W'Qji^sdWF-nc) 

LB(l+iaFR*[£b]Cm[0]) 

.  I  g*[4]\  1 g5T0p[4]ac^*[4]^(nf- 2^c) 


^na 


\^op[43/  '  ®0'  Lb{\  +  /ftF/^[£fc]CQW[0]){l  +  l(ft/?“ ^c)^*[^]^Qw[0]} 

Here,  l/Rd_£b]+  l/Rlp[£b]  with  l//?QP[^]  =  Pe[^]cr0p[w,  £b\  4>0CqW[0].  Moreover,  we  take  the  Fourier  trans¬ 

form  of  Eq.  (26)  and  get 


ftph(fiF)=ft°,{  S(dF)~  Ufif) 


From  Eqs.  (28)  and  (29),  we  find  the  drop  of  the  photoresponsivity  to  be 
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Aftph(oc)  =  k%- |_+J  <mF  nph(dF)  =^(L~^°])  I  f_+J  ddF  s^dF) 


0  | 


»2D  /  [  vi/(^*[4]cQW[o])z+n;j 

/[i-(/?*[4]/<[4])(^m/^o)]2+4(/?n4]cQW[03)2^ 

l  +  4(/?*[4]CQW[0])2n2 


which  is  zero  for  flc=0  and  n0ph(^mTo^£b]/n2D)  for  flc  III.  NUMERICAL  RESULTS  AND  COMPARISONS 
^l/(tf*[£*]CQW[0]).  The  chopping  frequency  at  which  the 

photoresponsivity  has  dropped  halfway  between  the  maxi-  ^  sec^on’  we  present  some  numerical  results  to 

mum  and  minimum  values  is  found  to  be  demonstrate  the  compensation  of  CF  in  quantum  wells  in  the 


ft*=l/(/?*[5,]CQW[0]).  (31) 

Figure  1  is  taken  from  an  early  measurement  performed 
by  Arrington,  et  al?  and  shows  the  quantum  efficiency-gain 
product  in  (a)  as  a  function  of  flc  and  the  extracted  in 
(b).  In  Fig.  1(a),  the  dynamical  suppression  of  the  photore¬ 
sponsivity  is  clearly  seen  as  a  function  of  £lc  for  various 
average  incident  photon  irradiances  <5  =  V([^Op(0]2)-  In 
Fig.  1(b),  is  plotted  as  a  function  of  <5  at  7=40  K  for 
different  dc  bias  voltages  Vb .  For  Vb-  - 1.5  V,  appears 
as  a  linear  function  of  <E.  When  Vb  is  increased  to  —3.5  V, 
we  find  is  a  linear  function  of  <I>  only  when  <I>  is  large 
but  becomes  independent  of  <5  when  <5  is  small.  We  know 
that  at  low  biases,  Rj_£b]>R°0V[£b]  and  1  fR*[£b~\ 
**  l//?op[£fc]  where  Rlp[£b]  is  proportional  to  <f>.  This  leads 
to  shown  in  Fig.  1(b)  as  predicted  by  Eq.  (31).  At 

high  biases,  the  situation  is  less  straightforward.  When  is 
small,  we  find  that  /?*[£*]  ^opC^]  and  is  independent 
of  <f.  At  higher  flux,  ^th  is  again  Hn- 

early  dependent  on  <5.  Consequently,  Eq.  (31)  is  verified 
experimentally,  which  proves  the  validity  of  Eqs.  (19)  and 
(20),  as  well  as  the  dynamical  photoresponsivity  7£ph(f)  de¬ 
fined  in  Eq.  (25). 

In  order  to  further  verify  the  relationship  in  Eq.  (32) 
beyond  the  limitation  of  <£m  /d>0<^  1 ,  we  present  in  Fig.  2  the 
numerical  results  for  A7£ph(flc)/7£ph  with  4>m/O0=l  and 
various  optical  fluxes,  device  temperatures  and  bias  voltages. 
From  the  solid  curve,  we  find  a  rolling-off  feature  of 
A7£ph(flc)/7£ph  at  ft^=0.3  Hz.  As  4>m  decreases  from  5 
X1012  cm"2  s"1  to  1X1012  cm"2  s"1  (dashed  curve),  the 
rolling-off  point  is  shifted  down  because  1  /R*[£b] 
~  l/ftop[£fc]  is  reduced  as  predicted  by  Eq.  (32).  When  Te 
increases  from  40  to  50  K  (dash-dotted  curve),  is  shifted 
up  because  l/R *[£b]  ~  VRJl^]  increases  with  Te  rapidly  as 
seen  from  Eq.  (32).  Finally,  when  we  reduce  £b  from  25  to  5 
kV/cm  (dash-dot-dotted  curve),  we  find  a  down-shift  of  fl  J 
since  l/R^'^h^l/R^h]  decreases  with  £b  as  implied  by 
Eq.  (32).  All  of  these  features  found  in  our  calculations  are 
consistent  with  the  experimental  observations  in  Figs.  1(a) 
and  1(b). 


I0W  „  10U 

Average  ImxUance  (ph/s-cm") 


FIG.  1.  Measured  quantum  efficiency-gain  product  of  a  sample  (see  Ref.  3) 
at  7^—40  K  and  central  wavelength  \0=7.1  fxm  as  a  function  of  the  chop¬ 
ping  frequency  in  (a)  for  different  average  photon  irradiances  <1>,  and  the 
extracted  chopping  frequency  at  which  the  quantum  efficiency-gain 
product  reaches  halfway  between  the  maximum  and  minimum  values  as  a 
function  of  4>  in  (b)  for  different  bias  voltages  Vb .  In  (a),  Vb  is  fixed  at 
—  2.5  V  and  $  is  varied  between  3.4X  1013  and  2.  IX  1015  ph/s  cm2.  In  (b), 
the  halfway  chopping  frequency  with  three  various  bias  voltages  Vb 
=  - 1.5,  —2.5,  and  —3.5  V  are  displayed  as  a  function  of  <$. 
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Frequency  &c  (Hz) 

FIG.  2.  Comparison  of  the  average  photoresponsivity  V{[^ph(nf-)]2)/7^  |Jh 
as  a  function  chopping  frequency  Clc  for  various  optical  fluxes,  device  tem¬ 
peratures  and  bias  voltages.  Here,  we  set  d>m  /<I>0=  1  in  our  calculation.  The 
solid  curve  corresponds  to  <J>m  =  5X  1012  cm-2  s“\  £b=25  kV/cm  and  Te 
=  40  K;  the  dashed  curve  to  1 X 1012  cm-2  s”1  but  the  same  Sb  and 
Te ;  the  dash-dotted  curve  to  Te-50  K  but  the  same  d>m  and  £b ;  and  the 
dash-dot-dotted  curve  to  £b  =  5  kV/cm  but  the  same  d>m  and  Te .  Other 
parameters  have  been  given  in  Tables  I-m. 

presence  of  both  £b(t)  and  4>op (t).  The  nonlinear  effect  be¬ 
yond  the  linear  model  will  be  considered  thereafter. 

In  our  numerical  calculation,  we  have  assigned  the  fol¬ 
lowing  forms  to  £b(t)  and  <l>op(0: 

eb{t) 

£0+Sm  sin[(277 t/Tp)  +  a0]  ,  jTp^t^(j+l/2)Tp 
0  ,  others  (32)  ’ 

for  7  =  0,  1,2,  ...  and 

3>oP(0 

O0+4>m  sin[27 rtlTp]  ,  jTp* (j  + 1/2) Tp 
0 ,  others  5 

(33) 

where  a0  is  the  phase  difference  between  £b(t)  and  <J>op(r) 
and  Tp  —  2irl£ic  is  the  time  period  for  both  £b(t)  and 
3>op(0.  £o  £m  are  the  dc  component  and  the  amplitude 
of  the  ac  component  of  £b(t).  4>0  and  <S>m  are  the  back¬ 
ground  flux  and  the  amplitude  of  the  modulation  component 
of  3>op(f).  Here,  we  take  rz=  1  ps.  The  other  sample  param¬ 
eters  used  in  our  calculation  are  listed  in  Tables  1-ID,  where 
the  photon  energy  is  in  resonance  with  the  energy  separation 
between  the  ground  and  excited  states  of  a  QW. 

A.  Compensation  of  charge  fluctuations 

In  this  part,  we  present  the  numerical  results  in  Figs.  3-5 
from  the  linear  model  in  Eq.  (24). 

TABLE  I.  List  of  internal  parameters  for  GaAs/ Al0  3Gao  7As  MQWs. 

Lw  Lb  «2D  &  vs  £s 

(A)  (A)  ( 10”  cm-2)  (10“4  cm-2)  (106  cm/s)  (kV/cm)  y  (meV) 

50  500  8  2.25  2  2  1 


# 


TABLE  H.  List  of  external  parameters. 


hco  (meV) 

£0  (kV/cm) 

a0  ( 1011  cm  2) 

T„  (s) 

Te  (K) 

168 

25 

7T  5 

400 

40 

Figure  3  presents  the  numerical  solutions  of  Eq.  (24)  for 
£na(0  as  a  function  of  time  t  for  different  values  of  £m ,  the 
ac  component  amplitude.  When  £m~0  (dash-dot-dotted 
curve),  due  to  CF  in  the  quantum  wells,  the  induced  nona- 
diabatic  space-charge  field  £na(t)  is  driven  by  d<&op(t)/dt 
over  the  first  half  period.  This  is  followed  by  a  decay  in  the 
second  half-period  when  $op(r)  — $0.  The  deviation  of 
£na(0  from  zero  measures  the  effect  of  the  CF,  which  is 
gradually  reduced  with  the  increase  of  £m  from  zero.  Be¬ 
cause  the  phase  difference  a0='ir,  £b(t)  and  <E>op(0  are 
completely  out-of-phase  with  each  other  and  the  effects  of 
the  CF  from  the  dual  tunneling  and  photon-assisted  escape 
are  compensated,  leading  to  a  decrease  of  the  deviation  of 
£na(0  from  zero. 

We  show  in  Fig.  4  the  total  nonadiabatic  photoemission 
currents  Ie[£b(t),  ®op(t)]-AIe(t)  as  a  function  of  time  t 
for  different  values  of  £m .  Within  the  adiabatic  limit  with 
£m~0  (dotted  curve),  the  photoemission  current  is  simply 
^[^o^opML  ^  it  is  proportional  to  <E>op(0*  Beyond  the 
adiabatic  limit,  the  total  nonadiabatic  photoemission  current 
is  dramatically  reduced  compared  to  the  adiabatic  value  in 
the  first  half  period,  but  is  enhanced  in  the  second  half  pe¬ 
riod.  This  is  due  to  the  fact  that  /^[^(O^opCO^A/^O 
always  remains  out-of-phase  with  d<£>op(t)/dt.  When  an  out- 
of-phase  electric  field  £b{t)  is  applied  to  the  system,  A Ie(t) 
which  results  from  CF  in  the  quantum  wells  in  both  the  first 
and  second  half  periods,  becomes  smaller  with  increasing 
£m .  We  attribute  this  to  the  cancellation  of  CF  in  the  quan¬ 
tum  wells  from  the  dual  tunneling  and  photon-assisted  es¬ 
cape  paths. 

The  results  of  the  dynamical  photoresponsivity  7£ph(0 
from  Eq.  (25)  as  a  function  of  t  for  different  values  of  £m  are 
compared  in  Fig.  5.  In  the  adiabatic  case  with  £m=0 ,  7£ph(0 
(dotted  line)  equals  its  static  value  7vph ,  and  is  independent 
of  3>op(0-  In  the  nonadiabatic  case,  lZph(t)  dramatically  de¬ 
creases  compared  to  7£ph  in  the  first  half  period,  but  in¬ 
creases  in  the  second  half  period.  In  the  first  period  where 
^na(0>0,  the  total  nonadiabatic  photoemission  current  is  re¬ 
duced  compared  to  its  adiabatic  value  as  shown  in  Fig.  4. 
However,  £na(t)  switches  its  sign  in  the  second  half  period  as 
shown  in  Fig.  3,  which  increases  the  photoemission  current. 
These  two  factors  together  give  rise  to  the  features  of  7Zph(t) 
seen  in  Fig.  5.  When  an  out-of-phase  electric  field  £b(t)  is 
applied  to  the  system,  7Zph(t)  gradually  approaches  71  ph  with 
the  increase  of  £m  due  to  the  compensation  of  the  CF.  We 
have  noted  that  the  dynamical  change  of  7Zph(t)  at  £b(t) 

TABLE  in.  List  of  emission  and  capture  parameters. 

SK  (kV/cm)  A0  £cp  (kV/cm)  Ba 

4.73  36  2.67  11.5 


FIG.  3.  Nonadiabatic  space-charge  field  EJf)  as  a  function  of  time  for 
different  ac-component  amplitudes  £m=0,  1,  2,  and  3  kV/cm. 

=  £0  reaches  approximately  as  high  as  40%  of  its  adiabatic 
value  7£ph.  However,  this  dynamical  change  is  greatly  sup¬ 
pressed  by  the  nonlinearity  as  discussed  next. 

B.  Effect  beyond  linear  model  . 

In  this  part,  we  present  the  numerical  results  in  Figs.  6-9 
which  compare  the  linear  model  in  Eq.  (20)  with  the  nonlin¬ 
ear  model  in  Eq.  (19)  at  £b(t)  =  £0. 

Figure  6  compares  the  calculated  total  electric  field  £0 
+  £na(f)  as  a  function  of  t  from  both  the  nonlinear  model  in 
Eq.  (19)  and  linear  model  in  Eq.  (20).  In  the  adiabatic  limit, 
the  total  electric  field  (dash-dot-dotted  line)  is  simply  £0.  The 
linear  model  introduces  a  large  CF  in  the  quantum  wells, 
which  enhances  £na(t)  in  the  system.  The  nonlinearity 
greatly  suppresses  the  CF  (solid  curve)  resulting  from 
d<bop(t)/dt*0,  giving  rise  to  only  a  small  deviation  of  the 
total  electric  field  from  £0.  We  know  that  the  positive  £aa(t) 
indicates  the  shift  down  of  the  Fermi  energy  E f  which  results 
from  charge  depletion  in  the  quantum  wells.  This  will  reduce 
the  dynamical  QW  capacitance  CQW[£na(r)]  compared  with 
CQW[0]  at  finite  temperatures  but  not  at  Te  =  0K  if 
From  this  figure,  we  see  that  the  dis¬ 
charging  process  [d£m(t)/dt>0]  is  accelerated ’in  the  non- 


FIG.  4.  Total  nonadiabatic  photoemission  current  Ie[£b{t),  <£op(0] 

—  AIe(t)  as  a  function  of  time  for  different  ac-component  amplitudes  £m 

—  0,  1,  2,  and  3  kV/cm.  For  convenience,  the  adiabatic  photoemission  cur- 
rent  Ie[£0l  <£<^(0]  with  £m=0  (dotted  curve)  is  also  included  in  the  figure. 


FIG.  5.  Dynamical  photoresponsivity  1 Zph(t)  as  a  function  of  time  for  dif¬ 
ferent  ac-component  amplitudes  £m=0,  1,  2,  and  3  kV/cm.  The  static  pho¬ 
toresponsivity  7£ph  with  £m=0  (dotted  curve)  is  also  included  in  the  figure 
for  the  sake  of  convenience. 


Unear  model  due  to  CQW[fna(f)]<CQW[0]  in  Eq.  (8),  which 
produces  a  small  discharging  constant  in  Eq.  (19). 

The  normaUzed  ratio  of  the  CF  in  quantum  weUs,  1 
~[AQ(t)/ £<S«2d]>  with  £b(t)~£ o  using  the  nonUnear  model 
is  presented  in  Fig.  7.  The  whole  fluctuation  process  within  a 
period  can  be  described  in  three  successive  steps  (see  Fig.  6): 
(1)  the  initial  discharging  process  d£na(t)/dt> 0,  (2)  the 
intermediate  charging  process  d£oa(t)/dt< 0,  and  (3)  the  fi¬ 
nal  discharging  process  d£aa(t)/dt> 0.  The  maximum 
“charge  depletion”  in  quantum  wells  is  about  10%  of  «2D  in 
the  initial  discharging  process  and  the  maximum  “charge 
accumulation”  is  15%  of  «2D  at  the  end  of  the  intermediate 
charging  process. 

Because  the  dynamical  photoresponsivity  TZph(t)  only 
depends  on  the  total  nonadiabatic  photoemission  current,  we 
show  in  Fig.  8  the  results  of  the  calculated  Ie[£b ,  <J>op(r)] 
-AIe(t)  as  a  function  of  t  with  £b(t)=£0  using  both  the 
nonUnear  model  in  Eq.  (19)  and  the  Unear  model  in  Eq.  (20). 
In  the  adiabatic  case,  the  photoemission  current 
(dash-dot-dotted  curve)  is  proportional  to 
®op(0*  Although  the  CF  in  the  quantum  wells  induces  a 


FIG.  6.  Total  nonadiabatic  electric  fields  £0  +  £na(t)  as  a  function  of  time 
from  the  nonlinear  model  in  Hq.  (19)  and  linear  model  in  Eq.  (20),  respec¬ 
tively.  The  dc  electric  field  £0  (dash-dot-dotted  line)  is  also  included  in  the 
figure  as  a  comparison. 
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FIG.  7.  Normalized  ratio  of  charge  fluctuation  [1-A QiOfeSn^]  calcu¬ 
lated  from  Eq.  (6)  using  the  nonlinear  model  in  Eq.  (19)  as  a  function  of 
time  for  Sb(t)=£0. 


FIG.  8.  Total  nonadiabatic  photoemission  current  Ie[£b(t),  ^op(r)] 
“A Ie(t)  from  the  nonlinear  and  linear  models  as  a  function  of  time 
for  £b(t)-£0.  The  adiabatic  photoemission  current  Ie[£b ,  <£>^(0]  with 
£b(t)=£0  (dash-dot-dotted  curve)  is  also  included  in  the  figure  for  the 
comparison. 


largevalue  of  A Ie(t)  in  the  linear  model,  it  is  dramatically 
suppressed  by  the  nonlinearity.  This  results  in  the  time- 
dependent  photoresponsivity  IZ^ h(t)  (solid  curve)  approach¬ 
ing  its  static  value  1Z ph  (dash-dot-dotted  line),  as  can  be  seen 
from  Fig.  9.  From  Fig.  9  we  also  see  that  the  dynamical 
change  of  TZ^(  t )  compared  with  7£ph  has  been  greatly  sup¬ 
pressed  by  the  nonlinearity  from  40%  to  10%.  Moreover,  the 
discharging  process  is  finished  with  a  much  higher  rate  due 
to  CQw[£na(r)]<CqW[0]. 


IV.  CONCLUSIONS  AND  REMARKS 

In  conclusion,  by  deriving  the  general  dynamical  equa¬ 
tion  for  Sna(t)  in  Eqs.  (23)  and  (24)  in  the  presence  of  both  a 
time-dependent  electric  field  £b(t)  and  a  time-dependent  in¬ 
cident  optical  flux  4)op(0  which  provides  dual  tunneling  and 
photon-assisted  escape  channels  to  the  system,  we  have 
found  a  compensation  of  the  CF  in  quantum  wells  when 
£b{t)  and  3>op(£)  become  out-of-phase  with  each  other.  By 
working  beyond  the  linear  model  in  Eq.  (20),  we  have  found 
from  a  nonlinear  model  in  Eq.  (19)  a  large  suppression  of  the 
nonadiabatic  deviation  of  the  photoresponsivity  and  a 
speed-up  of  the  discharging  process  due  to  the  depletion  of 
charge  in  the  quantum  wells.  For  a  special  case  with  a  dc 
bias  voltage  and  a  time-dependent  incident  optical  flux,  the 
experimentally  observed  drop  of  the  photoresponsivity  as  a 
function  of  the  chopping  frequency  has  been  reproduced  suc¬ 
cessfully  by  our  theory.  The  compensation  of  the  CF  in 
quantum  wells  predicted  in  this  article  will  be  verified  by  a 
future  experiment. 
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